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Abstract 

We study random-matrix ensembles with a non-Gaussian probability distri- 
bution P{H) ~ exp(— iVtr V{H)) where N is the dimension of the matrix H 
and V{H) is independent of N. Using Efetov's supersymmetry formalism, we 
show that in the limit N — > oo both energy level correlation functions and cor- 
relation functions of 5-matrix elements are independent of P{H) and hence 
universal on the scale of the local mean level spacing. This statement applies 
to each of the three generic ensembles (unitary, orthogonal, and symplectic). 
Universality is also found for correlation functions depending on some exter- 
nal parameter. Our results generalize previous work by Brezin and Zee [Nucl. 

Phys. B 402, 613 (1993)]. 
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The energy levels of a variety of physical systems including complex nuclei, disordered 
conductors and classically chaotic systems exhibit universal behaviour: The statistical prop- 
erties of the spectrum depend only on the fundamental symmetries of the underlying Hamil- 
tonian and can be described by random- matrix theory |1J. Three symmetry classes ex- 
ist: Systems with broken time-reversal symmetry are described by the unitary ensemble, 
time-reversal invariant systems by the symplectic or the orthogonal ensemble depending on 
whether spin-orbit coupling is present or not. The symmetry of the Hamiltonian does not 
specify the random-matrix ensemble completely. In general, one requires in addition that 
the matrix elements be statistically independent from each other. This additional condition 
restricts the probability density for all three ensembles to a Gaussian form and defines the 
Gaussian matrix ensembles. However, the assumption of statistical independence is used 
mainly for the sake of mathematical simplicity and is not motivated by first principles. 
Moreover, while accounting for the local energy level correlations on the scale of the mean 
level spacing d, the Gaussian ensembles fail to describe the global properties of the experi- 
mentally observed spectra. In particular, they predict a universal form for the mean density 
p{E) of eigenvalues (the semi-circle law for systems with a large number of levels) , whereas 
this quantity is known to be system specific and non- universal. 

Within random-matrix theory, more realistic forms of p{E) are obtained when one con- 
siders non-Gaussian ensembles. Based on numerical studies of such ensembles @ it has been 
conjectured |]] that their local spectral properties are independent of the measure, identi- 
cal to those of the Gaussian ensembles and hence universal. The effect of a non-Gaussian 
measure on the spectral statistics was first analytically investigated by Brezin and Zee [[3]. 
Within the unitary ensemble these authors proved universality for the local two-point level 
correlation function. (Properly smoothed level correlations on scales large compared to d 
were found to depend on the measure only through the endpoints of the spectrum. This 
latter result was later generalized to all three ensembles by Beenakker Q.) 

In the present paper, we prove the universality of arbitrary local correlation functions 
for any of the three generic matrix ensembles in the limit N — > oo, where N denotes the 
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number of levels. More specifically, we consider a correlation function C involving an arbi- 
trary number of level-density factors and/or S'-matrix elements. This function may depend 
parametrically on energy arguments and/or on additional parameters like the strength of 
an external magnetic field. We do allow for symmetry breaking. (An external magnetic 
field, for instance, breaks orthogonal symmetry). We compare the correlation function Cq 
evaluated for the Gaussian ensemble, and its analogue Cp evaluated for a non-Gaussian 
ensemble having the same symmetry and a distribution P{H) defined as in Eq. (D below. 
Under the single assumption that all parameters in C range over an interval A containing 
on average a finite number of levels A ~ 0(N~ l ), we show that Cg{<Ig{E)) = Cp(dp(E)). 
Here, dc{E) and dp(E) are the average local mean level spacings for the Gaussian and the 
non-Gaussian ensemble, respectively, both evaluated at energy E, which is located at the 
center of the interval A. These very general results are obtained by using the supersymme- 
try method. To allow for non-Gaussian probability measures, our derivation differs from the 
usual formulation of this method 

For definiteness we consider the unitary ensemble in the following. We emphasize that the 
orthogonal and the symplectic cases can be treated along exactly parallel lines. Accordingly, 
we study an ensemble of N x N Hermitian random matrices H with volume element d[H] = 
Yli>j dReHij ni>j dlrnHij. The probability density P(H) is defined by 

P(H) = Z- 1 exp{-NtrV(H)}, (1) 

where Z is a normalization constant. This is the most general density compatible with the 
basic assumption of random-matrix theory, namely that P(H)dH is invariant under unitary 
transformations H — > UHU^ 1 . The function V is assumed both to confine the spectrum to 
some finite interval and to generate a smooth mean level density, in the limit iV — > oo [0. 
Then, for the ensemble defined by Eq. ([I]), the mean level spacing d is of order iV -1 . Note, 
that V(H) = gH 2 , g > defines the Gausssian unitary ensemble. 

In the supersymmetry method, we generically express correlation functions as derivatives 
of a generating functional ||. The latter is written in terms of an integral over a supervector 
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\[> with bosonic (commuting) and fermionic (anticommuting) components 

I = J d[*](exp{^* t L 1 / 2 (H - E + M)L 1 / 2 \I>}), (2) 

where the brackets denote the ensemble average (...) = J d[H]P(H)(. . .). For example, 
in the case of the 2-point function at energies E and E' one defines the supervector by 
\& T = (Sf , Xi, $2, X2 ) with complex bosonic entries Si, S 2 and complex fermionic en- 
tries Xii X2, each entry being itself an ^-dimensional vector. The measure has the form 
d[ty] = Y[^=iltj=iidS*jdS l j_jdxl l jdXiJ,j- The Hamiltonian H is the direct product of the 
N x N Hamiltonian H and the unit matrix in the superspace. The energy E stands for 
the product of the mean energy E = (E\ + Ez)/2 and the unit matrix in both the level 
space and the superspace, L is the direct product of the unit matrix in level space with 
L = diag(l, 1, — 1, — 1). The matrix M contains energy differences and the source terms. 
In the case of scattering problems, M also contains couplings to external channels. To ac- 
count for dependences on external parameters, M may contain additional random matrices 
besides H over which additional ensemble averages must be performed. We postpone this 
calculation and confine attention to the ensemble average over H. In any case, we have 
M = 0(N~ 1 ) because we are interested in correlations involving energies of the order of the 
mean level spacing d ~ N^ 1 . For n-point functions with n > 2, the form of Eq. @ remains 
unchanged, although the dimensions of the vectors and of the matrices H, E, L, M 

in superspace will increase. Our proof applies to all these cases because it is independent of 
these dimensions. 

In the Gaussian case, one usually decouples the interaction generated by the ensemble 
average by means of a Hubbard- Stratonovich transformation ||. This procedure introduces a 
supermatrix a and maps the generating functional onto a non-linear a-model. The procedure 
relies on the Gaussian form of the probability density and does not apply to general P(H). 
The central point of our argument is based on the observation that it is nevertheless possible 
to introduce the "composite variables" a for any P(H). Indeed, for any P{H) the unitary 
invariance of the ensemble implies that for iV — > 00, the integrand in Eq. (g) depends on 



\P and ^ only via the invariant form A a p = j^L 1 ^ J2^=i ^ m^uS^/p ■ Here /i, v are level 
indices and a, (3, 7, 5 superindices ||. We explicitly introduce a supermatrix a with the 
same dimension and symmetry properties as A by writing / as an integral over a 5-function, 

1 = J d[V] J dad(a - A)(exp{ i -^L 1/2 GL 1/2 ^}}, (3) 

with the abbreviation G = H — E + M. The (^-function is replaced by its Fourier represen- 
tation 

I = J d[^f] J da J dr exp{^iVtrg(Yo-)} 

x(exp{^ t L 1 / 2 (G-r)L 1 / 2 ^}), (4) 

and the multiple Gaussian integral over the \l/-supervector is performed 

1 = JdaJdTexp{^Ntig(r(r)}(exp{-^tTtTg ln[G - r]}). (5) 

We have now expressed the functional J as a integral over two coupled supermatrices a and r 
which contain all relevant degrees of freedom. In the limit iV — > 00, the remaining integrals 
can be done explicitely using the saddle-point approximation. In particular, it will turn out 
that the (diagonal) saddle-point of the cx-integral determines the mean level density. These 
steps will prove the claimed universality by comparison with the well-known Gaussian case. 

To perform the ensemble average we transform H to diagonal form, (UHU^) = A, and 
integrate separately over eigenvalues A and eigenvectors U. Expanding in powers of M, we 
have 

(exp{ — -trtrg ln[G — r]}) = (exp{ — -trtrg \\\D} 

xexp{ — trtrg [1 + ]T ^ '- (D- 1 WMU)) n ]}), (6) 

~~ n=l ^ 

where D = (A — E — r) is diagonal in the level space. The expansion in powers of M in 
Eq. (^|) cannot be terminated with the first-order term because any power of M may be 
of the same order A^ _1 as M. This is the case for S'-matrix correlation functions ||. The 
distribution of eigenvectors does not depend on the form of the probability density in Eq. 



(H). In the large- N limit the eigenvectors are Gaussian distributed and the average over 
eigenvectors is evaluated using Wick contractions. To leading order in powers of N' 1 we 
find that the last exponential in Eq. (|6|) takes the simple form — |trtrg[l + (-^tr.D _1 )M]. 
The remaining eigenvalue-integrations are done using the saddle-point approximation |jTD| . 
Explicitly, the average over eigenvalues appearing in Eq. (|6|) involves the exponential 

- \ £ trg In (£>)„ - ~tr trg In 1 + ( 1 £(ZT ^ M 

-i\r£V(A M ) + 2£ln|A M -A,|, (7) 

where the last term results from the Jacobian associated with the transformation from the 
matrix elements of H to its eigenvalues. Of the four terms in expression (0) the first is 0(N) 
and the second 0(1). The last two terms are 0(N 2 ) and determine the saddle point values 
A^ p . The calculation is explicitly carried out in ref. |TTJ and introduces the average local 
level density p(E), and the resolvent F, defined by 

h E ~ / dE '^T^ = F(E + r). (8, 

Substituting A^ p for A M in the first two terms of Eq. (0), expanding the terms ~ 0(N 2 ) 
around the saddle-point and performing the Gaussian integrals, we find 

/ =JdaJdr exp | ^iVtrg(ra) -^trg ln[A^ p - E— r] j 

x exp j— - tr trg ln[l + F(E + r)M]| . (9) 

(The integration over the A M around the saddle-point cancels against the normalization Z.) 
Again, the saddle-point approximation is used to integrate over r. The last exponential in 
Eq. (^D has a term of order 0(1) in the exponent and can be omitted. Hence, for fixed 
a the equation ia = F(E + r sp ) determines the saddle-point r sp (cr, E). By expanding the 
exponent to quadratic order in the fluctuations St, one can easily verify that the integral 
over St yields unity. Therefore, we obtain 

J^expj^iVtrg^) - ±£trgln(A*-£-7 J *)} 

x exp|--trtrgln(l - iaM)\ , (10) 



where now only the integration over the supermatix a remains to be done. The saddle-point 
a sp is found from the first two terms in the exponent 

<9r sp <9r sp 
^ sp ^U + ir s * = F(E + ^ SP )^U • (11) 

The saddle-point equations for r and a together show that r sp (a sp ) = and ia sp = F(E). 
One observes that the saddle-point equation for a sp is invariant under pseudo-unitary trans- 
formations. This implies that also <jq = T _1 a sp T is a saddle-point, where T generates 
pseudo-unitary transformation on the space of supermatrices. In general <tg 7^ ^ and 
hence the ctq form a manifold of solutions of the saddle-point equation. We expand the 
exponential in Eq. (|T0|) in the vicinity of the saddle-point solution ctq, carry out the integral 
over the massive modes (which gives unity) and find the result 

I = J dfi(t) exp {--trtrg ln[l - T~ 1 F(E)TM\\ , (12) 

where the integration is now over the manifold of saddle-points. As usual one has to 
give E an imaginary part such that ImF(E) ~ L to guarantee convergence. Both the 
structure of the saddle-point manifold and the measure dp depend only on the symmetry of 
the ensemble and on the dimension of the supervectors W in Eq. (^). In particular, both 
are independent of the probability density P{H). On the other hand, the latter specifies the 
mean level density p(E) and the function F(E) and hence sets the local energy scale. To 
see this most clearly one chooses E such that KeF(E) = and therefore F(E) = —iirp(E)L 
0. For a symmetric confining potential V one may take E at the center of the spectrum, 
£7 = 0. Then 

I = J d/i(t)exp j--trtrg ln[l + mp{E)T- 1 LTM] \ . (13) 

This proves universality: Our result has the same form as in the Gaussian case 0, the 
probability density enters only through the local mean level density p(E). On the scale 
set by p{E) all derivatives of / and hence all correlation functions are independent of the 
probability density P(H) and are thus universal. 



The proof, presented here for the unitary ensemble, applies equally to the orthogonal 
and the symplectic ensemble. In either of the latter cases, the structure of the supervector 
\l/ differs from the unitary case. This structure is reflected in the symmetry properties 
and dimensions of the matrices a and r and, eventually, of the matrix T which generates 
transformations on the saddle-point manifold. However, our threefold use of the saddle-point 
approximation is completly independent of such symmetry properties. 

Last we turn to correlation functions which depend on an external parameter ||11|| . One 
has to distinguish between two cases: The external perturbation may either preserve (case 
(i)) or violate (case (ii)) the symmetry of the original matrix ensemble. We consider class 
(i) first and demonstrate our point for the case of the 2-point level correlation function. It 



is defined by (tr(£ + e/2 - H - Ja/NH')tr{E - e/2 - H + yJa/NH')) where the energy 
difference e is 0(N~ l ) and H' is a matrix ensemble with the same symmetry and the same 
distribution function as the original ensemble H. The average is over the distributions 
of both H and H'. For simplicity we restrict ourselves to symmetric distributions. The 



factor y/N 1 appearing in the definition of the 2-point correlation function ensures that the 
corralations decay on the typical scale a ~ 0(1). In the generating functional of Eq. (0) 



there now appears an additional term ^a/NH'L which can be included into the definition 
of M by replacing M — > M + J~^H'L. After averaging over the ensemble H one finds that 
/ is given by 

/ = y*d//(t) (expj-^trtrg In [l + m p{E)T~ 1 LTlsJi\\) H > ■ (14) 



Now the logarithm is expanded in powers of H' . Due to the factor ya/iV, taking the ensem- 
ble average over H' reduces to calculating the second cumulant, all higher-order cumulants 
being small in comparison by at least a factor iV -1 / 2 . The final result for I is 

1 = J d//(t)exp j--trtrg In [l + m p(E)T' 1 LTM^ j 

x exp |-^(^) 2 [trgCT-^TL^tr (H')%] j. (15) 

The distribution over H' enters only trough its second moment. Its value for a non-Gaussian 



distribution for H' differs from what one would find for a Gaussian distribution. However, 
this difference does not affect the form of the correlation function and only leads to a 
rescaling of the parameter a. Case (ii) is treated along exactly parallel lines and leads 
to exactly the same conclusion: Aside from a scaling factor affecting the parameter which 
governs symmetry breaking, the form of the correlation function is the same for Gaussian 
and non- Gaussian ensembles. 

In summary, we have investigated the consequences of non- Gaussian probability measures 
within random-matrix theory. We have shown that in the limit of a large number of levels 
global and local properties of the spectrum separate. Global properties like the mean level 
density do depend on the form of the measure. Local properties, in contrast, are independent 
of the measure. They are determined only by the symmetry of the ensemble and they 
are identical to those of the corresponding Gaussian ensemble. This holds for all three 
generic ensembles and for arbitrary form of the measure. Our analytical result establishes 
generally and for the first time that all local random-matrix correlations are independent of 
the measure and hence universal. 

We thank F. von Oppen and J. Zuk for helpful and informative discussions. 
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